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Abstract —Energy harvesting (EH) cognitive relays are an exciting 
solution to the problem of inefficient use of spectrum while achieving 
green communications and spatial diversity. In a spectrum sharing 
scenario, we investigate the performance of a cognitive relay net¬ 
work, where a secondary source communicates with its destination 
over Nakagami-m channels via decode-and-forward EH relays while 
maintaining the outage probability of the primary user below a 
predefined threshold. Specifically, we derive a closed-form expression 
for the secondary outage probability and show that it is a function of 
the probability of an EH relay having sufficient energy for relaying, 
which in turn, depends on the energy harvesting and consumption 
rates of the EH relay and the primary outage probability threshold. 
We also show that relaxing the primary outage constraint may not 
always benefit the cognitive EH relay network due to the limitations 
imposed on the relay’s transmit power by the energy constraint. 

1. Introduction 

With the ever-increasing demand for wireless services along 
with a need for green communications, spectral efficiency and 
energy efficiency have become important criteria in the design of 
future wireless systems. Energy harvesting (EH) cognitive radio 
m-nii is a promising solution to improve the spectrum utiliza¬ 
tion; in particular, spectral efficiency is improved by spectrum 
sharing, while achieving self-sustaining green communications. 
In cognitive radio, a secondary user (SU) may share the spectrum 
with a primary user (PU) provided that the interference from it 
to PU remains below a given threshold El. 

The use of cooperative relays in cognitive radio has gained 
significant attention as they have the potential to improve the 
coverage and reliability of SU’s transmission while sharing the 
spectrum with PU |[6l- |[T6l . However, the relays may have 
limited battery reserves, and recharging or replacing the battery 
frequently may be inconvenient. This invokes the need for an 
external power source to keep relays active in the network. The 
EH relays can overcome such energy shortage while exploiting 
the spatial diversity IfTTl-ini. As to the EH relays in cognitive 
radio, Ea, ED consider an EH secondary relay which helps 
relaying the secondary data, and perform the secondary outage 
analysis for Rayleigh fading channels under the interference 
constraint at the primary receiver; while in ||24l, cooperative 
communication via multiple EH relays is considered. 

In this paper, we consider the case where SU uses the best relay 
from multiple EH relays for its own transmission over Nakagami- 
m channels, given that PU’s outage probability remains below a 
given threshold—we characterize the interference to PU by its 
outage probability. Eor EH relays, the optimal use of available 
energy is crucial. Low transmission power to conserve energy 
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may prolong the lifetime of a relay, however, at the cost of 
increased outage; whereas higher transmission power improves 
the transmission quality, but at the expense of higher energy 
consumption rate reducing the future chances of transmission. 
Due to this EH nature of relays, the best relay selection becomes 
tricky as only relays having sufficient energy to forward the data 
to the destination, called active relays, can be considered for the 
selection, making energy a crucial factor in the relay selection. 
Additionally, in spectrum sharing, the secondary communication 
via EH relays differs from that in non-spectrum sharing environ¬ 
ment; because, in spectrum sharing, EH relay’s transmit power 
depends not only on the energy availability with it, but also on the 
maximum power allowed by PU’s outage constraint. Eor example, 
in a case where a relay has harvested less energy, it may not trans¬ 
mit with the maximum power allowed by PU’s outage constraint. 
On the contrary, even if the relay has harvested large amount 
of energy, a tight PU’s outage constraint may not allow relay 
to transmit with higher power. Thus, there exists an interesting 
tussle between these two constraints putting a stronger restriction 
on relay’s transmit power than it would have been in the case 
with only one constraint, i.e., spectrum sharing without energy 
harvesting or non-spectrum sharing energy harvesting. Intrigued 
by the aforementioned tussle, in this work, we investigate its 
impact on the secondary network’s performance, which is missing 
in (m-jm. The main contributions of this paper are as follows: 

• Eirstly, with the best relay selection scheme that maximizes 
SU’s end-to-end signal-to-interference-noise-ratio (SINR), 
we specifically derive a closed-form expression for the 
outage probability of EH decode-and-forward (DE) cognitive 
relay network in Nakagami-m channels under PU’s outage 
constraint. We also consider the interference from PU while 
deriving the outage probability expression. 

• Secondly, for better utilization of the harvested energy, we 
calculate the probability of a relay being active. We show 
that, besides energy harvesting and consumption rates, the 
probability of a relay being active depends on PU’s outage 
probability threshold. We then couple the energy constraint 
due to the EH nature of relays with PU’s outage constraint. 
We investigate which of the two constraints dominates the 
performance of EH relays and find the respective regions of 
dominance that regulate the transmit powers of EH relays. 

• Einally, we investigate the effects of fading severity parame¬ 
ter, number of relays, and the average energy harvesting rate 
on the secondary outage probability and tradeoff between the 
energy constraint and PU’s outage constraint. 

H. Energy Harvesting and Spectrum Sharing Model 

As shown in Eig.[Tl the network consists of a primary trans¬ 
mitter (PT), a primary destination (PD), a secondary transmitter 



Fig. 1. Secondary transmissions via EH relays with underlay spectrum sharing. 

(ST), a secondary destination (SD), and M energy harvesting 
DF secondary relays (SRs). The PT, PD, ST, and SD are 
conventional nodes with constant energy supply (e.g., battery). 
The ST-SD direct link is assumed to be unavailable O-llIll- 
The ST communicates with SD over ith half-duplex EH relay 
(SRi), i G The channel between a transmitter 

p e {PT, ST, SR^} and a receiver q G {PD,SD,SR^}, is a 
Nakagami-m fading channel; hp-q denotes the channel coeffi¬ 
cient. Thus, the channel power gain \ hp-q\‘^ is Gamma-distributed 
with mean ftp-q and fading severity parameter rrip-q. We can 
write the probability density function and cumulative distribution 
function (CDF) of a random variable U = \hp-q\‘^ as 

TTip — q 

a ^ ^ 

fu(u) = (1) 

i [rrip-q) 
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respectively, where r('), r(', ^nd T(',') are the complete, 
upper incomplete, and lower incomplete Gamma functions 1^ . 
respectively; ap-q = mp-qjVtp_q. The channels are independent 
of each other. For PT-PD, ST-PD, and SR^-PD links, we assume 
the mean channel power gain knowledge due to limited feedback; 
while SRi and SD have the instantaneous channel gain knowledge 
for respective receiving links, i.e., ST-SR^ and PT-SR^ links 
at SR^, and SR^-SD and PT-SD links at SD (91, (261 . The 
secondary communication happens over two phases, each of 
T-second duration. All channels experience block-fading and 
remain constant for 2T-second, i.e., two phases of secondary 
communication, as in ini-ci. In phase 1, ST transmits to EH 
secondary relays, while in phase 2, the received signal from ST 
is forwarded by one of the relays to SD after decoding. Note that 
in phase 2, no relay might be active due to the lack of energy 
required to forward the signal. 

A. Energy Harvesting Model 

The energy harvesting process of a relay i is stationary and 
ergodic 113, with mean Ffav,i Joules per second (J/s). This 
model encompasses different energy harvesting sources like solar, 
vibrations, radio frequency (RE) in the surroundings (771 . and 
different energy harvesting profiles (28l. The EH relay stores 
the harvested energy in a battery with negligible leakage. For 
analytical tractability, we assume the capacity of the energy 
storage to be large ca -(20l . In addition, the energy consumption 
occurs only in data transmission; any other energy expenditure, 
e.g., energy consumption in signal reception and processing, is 
not considered for the purpose of exposition (TSl, (191 , (22l. 


B. Maximum Secondary Transmit Powers in Spectrum Sharing 

In this work, we characterize the quality of service (QoS) 
of PU by its outage probability. For constant transmit power 
of PT (Ppx), the PU outage probability should be below a 
certain threshold 0p given the interference from the secondary 
transmitter and the relay. This constraint limits the transmit 
powers of ST and SR to Pst and Psr, respectively. In phase 1, 
the outage probability of PU Pp,out,ST when ST is transmitting, 
is given as 


Pp,out,ST = Pr (log 2 (1 + 7pd) < Pp) < Op, 


( 3 ) 


where 7 pd = Pst7t-pdI°+^o 

noise power of additive white Gaussian noise (AWGN) at all 
receivers, and Pp is the desired data rate on the primary link. 
Proposition 1. We write Pp,out,ST as follows: 


p,out,ST 


= 1 - 


exp 


— apT-PD^p-Vp 

. ^PT 


mpT-PD — 1 

X E 

k=0 


r(mST-PD 

OpT-PD^pA^O 

PpT 

r(msT-PD + t) 


1 ^ 
hE 


k\ 


(t)' 


/^Q!PT-PD^pr"ST I PD+^ 

( -- + OST-PD ) 


(4) 


Proof: The proof is given in Appendix |Al ■ 

Using the value of maximum ST power Pst allowed by 
PU’s outage constraint 0p can be numerically found. Similarly, 
in phase 2, the maximum allowable transmit power Psr for a 
relay can be numerically found from (5]), replacing the role of 
secondary transmitter in (H) by the secondary relay and replacing 
corresponding channel parameters. 
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C. Active Relays and Best Relay Selection 

Assume that out of the total M relays, N relays are active 
{N G {0,1,...,M}) due to energy availability, and a relay has 
to be selected from active N relays. An active relay is the relay 
having sufficient energy to forward the received data from ST. 

For an opportunistic DF relaying, the relay with the largest end- 
to-end SINR at SD, called the best relay, is selected to forward 
the signal. When N relays are available for selection, the largest 
end-to-end SINR at SD is given by 

7t7 = max (min(7sR,, 7R; D)), (6) 

where M is the set of active relays. Note that M is an empty set 
when no relay is active. 7 sr. and 7r-d are SINRs at ith relay 
and at SD over R^ — D channel, respectively, and are given as 
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where we obtain Pst and Psr. from dH) and ([5]), respectively. 
III. Secondary Outage Analysis 

Now, when we select the best relay out of N active relays, the 
secondary outage probability P^out t)e given as 

P7ut(7) = Pr(7t7 < 7) = Pr(lmax^(min(7sR,,7RiD))<7Y 

1 . 

where secondary’s desired secondary rate is Pg = ^ lc)g 2 (1 + 7 )- 
For the ease of representation and without compromising 
the insight into analysis, we consider msT-SRi = R^st-sr, 
ttisr.-sb = msR-SD, rnpT-SRi — R^pt-sr, R^sr—pd = 
RT'SR-PD, ^ST-SRi =^ST-SR, ^SR^-SD = f^SR-SD, f^PT-SR^ = 

^PT-SR, and f^sRi-PD =^sr-pd- Then, Psr. = Psr. Below 
we give the closed-form expression for P^out- 

Proposition 2. We write P^o^t (7) follows: 
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where A is given by (1221) . 

Proof: See Appendix |Bl A key idea in the proof is to 
consider the dependency between 7 r.d and 7 r^d {i 7 ^ jA ^ 
{1,2 ,..., M}), originating due to the common term |/ipT-SDp- 


For an EH relay, its operation is subject to the energy neutrality 
constraint, which states that a relay cannot spend more energy 
than it has harvested. Thus, it is possible that the relay might 
remain inactive for some time due to the lack of energy. 

Let us denote the probability of a relay i being active by 7 ^ > 0. 
In a non-spectrum sharing scenario, r]i depends on the relay’s 
energy harvesting and consumption rates. Based on these two 
factors, the relay i operates in two regions as follows: 

• Energy constrained region (rn <1) 

• Energy unconstrained region (rji = 1). 

A relay operates in the energy unconstrained region if its av¬ 
erage energy consumption rate is less than the average energy 
harvesting rate, i.e., the relay is always active. 

We assume Pav,i = ^av without loss of generality. Then, we 
have r]i = 7]. The energy available with the relay depends on the 
factors that, how frequently the relay is selected; its harvested 


energy till now; and when was the energy harvested. As we 
will show later, in the case of spectrum sharing with PU, the 
probability of a relay being active, i.e., 7 , depends not only on the 
energy harvesting rate, the total number of relays in the system, 
the energy consumed by a relay in its each transmission, but also 
on PU’s outage constraint. Using the following proposition given 
in m, we show the dependency of 7 on PU’s outage constraint. 
Proposition 3. Let the probability of selecting a relay be uj. Then, 
uj = The relays remain active with the probability 

7 ^= 1 - [( 1 -Mcj)+]^. ( 11 ) 

All the relays become energy unconstrained, i.e., r] = 1, when 
UJ > 1/M. We denote (t)+ = max(0,T). 

The expression for uj in Proposition 3 is obtained from the 
energy neutrality constraint and stationarity and ergodicity of the 
energy harvesting process. From Proposition [3l one can notice 
that the probability of a relay being active depends on the power 
Psr with which the relay performs a transmission. Equations ^ 
and (HB together show that Psr, in turn, the probability of a 
relay being active, depends on the primary outage constraint. 

Given N out of M relays are active, each with the probability 
7 , we obtain the final expression for the secondary outage 
probability with EH relays by unconditioning over the number 
of active relays as 

]V[ / -J, 

p.„. = E ( 12 ) 

N=1 ^ ^ 

where P^out given by (fTOl) is the secondary outage probability 
when we select the best relay among active N relays to forward 
the signal from ST; Pgout i^ secondary outage probability 
when no relay is active, and is equal to 1 . 

IV. Discussions and Results 

In spectrum sharing, PU’s outage constraint governs the max¬ 
imum transmit power of relays. In addition, if relays are energy 
harvesting, due to the limited available harvested energy, the 
probability of a relay being active plays an important role in 
the performance of the secondary system. In this section, we 
will first discuss the effect of PU’s outage constraint on SU’s 
outage performance for the case when an EH relay on selection, 
uses the maximum power allowed by the primary, aiming to 
reduce the secondary outage probability; even though, it might 
also reduce the relay’s probability of being active. Then, we 
will consider the case when, along with PU’s outage constraint, 
EH relays aim to keep their probability of being active to one 
(7 = 1)—which we will call the energy constrain)^ In this case, 
we will show that the transmit power of the relay is regulated by 
the dominant of the two constraints—PU’s outage constraint and 
energy constraint—and we will find the region of dominance for 
each constraint. Finally, we will see from results that in both the 
above cases, relaxing PU’s outage constraint beyond a level does 
not offer any benefit to the secondary system with EH relays. 

^Note that the constraint 77 = 1 is different from the energy neutrality 
constraint. With latter, the consumed energy by a relay cannot exceed its harvested 
energy, whereas the constraint of always being active does not allow relay’s energy 
consumption rate to exceed its energy harvesting rate. With the higher energy 
consumption rate, the relays will eventually consume the energy harvested in the 
past before acquiring sufficient newly harvested energy to keep them active. 














A. System Parameters and Simulation Setup 

We consider following parameter values: PU transmit power, 
PpT = 15 dB; the desired primary rate, IZ^ = 0.4 bits/s/Hz; the 
desired secondary rate, IZ^ = 0.2 bits/s/Hz; noise power, Nq = 
—60 dBm. Denote fading severity parameters on forward and 
interference channels by mf and mint, respectively. We consider 
a 2-D topology, where {xi, yi) defines the coordinate of ith user. 
The mean channel gain between ith user with coordinate {xi, 
yi) and jth user with coordinate {xj, yj) is where dij is 

the distance between users i and j in meters and A is the path- 
loss coefficient which is assumed to be 4. Without any loss of 
generality, ST is placed at (0, 0), M relays are clustered and 
collocated at (50, 0), and SD is placed at (100, 0). Also, PT and 
PD are located at (50, 50) and (100, 50), respectively. 

B. Effect of PU's Outage Constraint 

Figs.[2|and[3]show the effect of PU’s outage constraint on the 
outage probability of the secondary system with EH relays 0 The 
selected EH relay transmits with the maximum power allowed 
by PU’s outage constraint. We notice that the increase in the 
primary outage threshold 0 p increases the maximum transmit 
power PsR allowed for the relay, which initially reduces the 
secondary outage probability Ps,out- However, with the increase 
in the threshold 0p beyond a level, a tipping point will be 
reached after which Ps,out will increase even with the increase in 
PsR as relays will consume energy at a higher rate than they 
will harvest, i.e., relays will become energy constrained (see 
plots for Pav = l 5 2 in Fig.©. This will reduce the probability 
of a relay being active, thereby reducing the number of relays 
available to forward the data to SD. As long as Psr is below 
a certain level so that uj = > 1/M as shown in Proposition 

3, the relays operate in the energy unconstrained region, i.e., the 
harvested power is more than the transmit power Psr. But, with 
relaxation of PU’s outage constraint, eventually the value of Psr 
increases such that uj < 1/M, making relays energy constrained 
and increasing Ps,out- Also, increase in the harvesting rate Pav 
delays the occurrence of the tipping point as expected, and at high 
harvesting rates, the relays might operate completely in the energy 
unconstrained region due to availability of abundant energy (see 
the plot for Pav = 4 in Fig.[2]). 

Remark 1. Relaxing PU’s outage constraint may not always 
improve the performance of EH secondary relays in spectrum 
sharing. That is, unlike conventional non-EH case, due to lack of 
energy, the relays with EH capability may not transmit with the 
maximum allowed power even though they are allowed to do so. 

C. Effect of Eading Severity Parameter 

Fig. [3] shows the effect of fading severity parameter on Ps,out- 
We notice that, before the tipping point, i.e, in the energy 
unconstrained region, Ps,out is lower for higher fading severity 
parameter mf on forward channels. However, the trend reverses 
after the tipping point. This is because, with the increase in mf, 
the fading effect subsides over the primary link between PT and 
PD, providing an extra margin for maximum secondary relay 

^Simulation results validate the analysis. The number of iterations is up to 10®. 



Fig. 2. Secondary outage probability (Ps,out) vs. primary outage probability 
threshold (0p), M = 3, mf = 2, mint = 1- 



Fig. 3. Secondary outage probability (Ps,out) vs. primary outage probability 
threshold (0p), effect of fading severity parameter and number of relays M, 
'^int — 1? -ffav — 2 J/s. 

transmit power Psr for a given 0p. This helps in achieving lower 
Ps,out for higher mf in the energy unconstrained region, where 
the energy harvesting rate is higher than the energy consumption 
rate. As shown in Fig.O for a given harvesting rate, due to higher 
allowed Psr (higher energy consumption rate) for higher mf, 
the tipping point arrives earlier than that for lower mf. After the 
tipping point, since relays enter the energy constrained region, 
higher mf, in turn, higher energy consumption rate, reduces the 
probability of a relay being active. This often leads to non¬ 
availability of relays for transmission, increasing Ps,out fof higher 
mf. Also, increase in the number of relays M increases the 
probability of being active (see CB) as the candidate relays for 
cooperation increases (increased diversity) due to which a certain 
relay is chosen less frequently. This reduces Ps,out- 

D. Joint Effect of PU’s Outage Constraint and Energy Constraint 

From discussions of Figs. [2| and [51 we note that EH relays being 
inconsiderate towards their probability of being active, leads to 
SU’s inferior outage performance beyond the tipping point. Now, 
for instance, assume that EH relays try to remain always active 
(j] = 1 ), i.e., try to satisfy the energy constraint, irrespective 
of PU’s outage constraint, and transmit with power Pg^a- From 
Proposition 3, we can see that satisfying the energy constraint 
corresponds to uj>l/M, i.e., Ps,a<2MPav That is, as long 
as EH relays transmit with power no greater than 2MPav, they 



































Fig. 4. Secondary outage probability (Ps,out) vs. primary outage probability 
threshold (0p), M = 3, = 3, mint = 2, ifav = 2 J/s. 

always remain active. Now, if we combine the energy constraint 
with PU’s outage constraint, Fig.[4| shows that in the energy 
unconstrained region, though an EH relay may transmit with 
maximum power 2MHav maintaining = 1, the power 2MHav 
does not satisfy PU’s outage constraint, i.e., 2MHav > Psr- 
This leads to higher Ps,out in EH relay spectrum sharing scenario 
governed by both the energy constraint and PU’s constraint than it 
would have been in EH non-spectrum sharing scenario governed 
by the energy constraint alone. In the energy constrained region, 
the energy constraint becomes dominant, i.e., 2MPav < Psr- 
Thus, even though PU’s outage constraint is satisfied, and allows 
EH relays to transmit with the maximum power Psr, the energy 
constraint is violated, causing r] < 1 and increasing in Ps,out as 
discussed for Eigs.[2|and[3l Therefore, we can see from the above 
discussion that, to satisfy both constraints, the maximum power 
with which EH relays may transmit is min(2MPav5 Psr)- As 
shown in Fig.[4l in the energy constrained region (after the tipping 
point), though transmitting with power min(2MPav5 Psr) avoids 
the increase in Ps,out, relaxing PU’s outage constraint does not 
improve SU’s outage performance. 

V. Conclusions 

Under the primary outage constraint, this paper has analyzed 
the outage performance of the secondary communication via 
energy harvesting relays in Nakagami-m channel. In a spectrum 
sharing scenario, the results show that, besides energy harvesting 
nature of relays, the primary outage constraint also strongly 
infiuences the probability of a relay being active. We note that 
relays should keep their probability of being active to one; 
otherwise, obeying only the primary outage constraint may lead to 
the inferior secondary outage performance. That is, in energy har¬ 
vesting spectrum sharing, due to the energy constraint, relaxing 
the primary outage constraint may not always improve secondary 
outage performance unlike in non-energy harvesting case. Further, 
we have found the region of dominance for each of the constraints 
and proposed the optimal transmit power for relays to subside 
their inferior performance in the energy constrained region. We 
observe that increase in the number relays and lower fading 
severity parameter delays the entry of relays into the energy 
constrained region, boosting the secondary outage performance. 


Appendix A 
Proof of 0 

From ([3]), conditioned on |/ist-pdP = x we can write 
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Then, using (O in ([T3]) and unconditioning over |/isT-PDp, we 
can write @ as 
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Simplifying ([T5l) and using binomial expansion, we get 
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Solving (O, we get the required expression in 0. 

Appendix B 
Proof of (fTol) 

From®, we can see that 7 r.d and 7 r^d 0 

{1,2,...,M}) contain the common term |/ipT-SDp, which 
makes them dependent. Thus, conditioned on \hpT- SDp = X, 
we can write the CDF of as 
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Using the procedure to derive Q, we can write (l2Tll as 
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Substituting (l20l) and (l22l) in (fTTI) and using the multinominal 
theorem ll25]| . we get 
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k{rk-i — Tk). In the following 
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Step, we use binomial expansion of (l + j 
take expectation over |/ipt-sdP- Then, we can write 
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(mpT-SD — 1 )! 

Solving the integration in (l24l) . we get the required closed-form 
expression for the secondary outage probability given by (Hii. 
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